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^ ■ Abstract 

, Problems in realization of silicon-based solid-state NMR quantum computer with ensemble 

addressing to qubits are considered. It is presented the extension of Kane's scheme to ensemble 
OO , approach version with strip gates. For the initialization of nuclear quantum states it is proposed 

to use the solid-state effect in ENDOR technique whereby the nuclear spins can be practically 
fully polarized or, that is the same, indirect cooled to the spin temperature less than ~ 1 mK. 
^ ' It is suggested the possible planar silicon topology of such ensemble quantum computer and 

, shown that the measurement with standard NMR methods signal of L ^ 10^ qubit system may 

OO ' be achieved for a number of ensemble components N > 10^. As another variant of ensemble 

silicon quantum computer the gateless architecture of cellular-automaton is also considered. The 
decoherence of quantum states in the ensemble quantum computers and ways of its suppression 
is also discussed. 
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Introduction 

^ . Atomic nuclei with a spin quantum number 1=1/2 seem to be the natural candidates for qubits 

^ : 



> 

^ : 

, — two- level quantum elements in quantum computers. The early approach to NMR quantum 
^ ' computers was suggested in 1997 [||, ^ and then confirmed in experiments [§, ^. In this approach 
it was used several diamagnetic organic liquids which individual molecules, having a small number 
of interacted non-equivalent nuclear spins-qubits with 1 = 1/2 and being nearly independent one 
^ ' another. They act in parallel as an ensemble of almost independent quantum microcomputers. 
■ Since in liquids the nuclear spins are very weakly coupled with the environment, the consider- 

ation may be restricted to nuclear spins of individual molecules {reduced quantum ensemble). For 
finite temperatures, the reduced ensemble is in the mixed state described by a density matrix. 

Initialization of the nuclear spin states in this case means the separation of blocks from density 
matrix that must have properties similar to those of pure states (a pure-state diagonal matrix 
can have only one nonzero element). The states described by such matrix are called effective or 
pseudo-pure states. Several way of initialization have been reported in [||, ^, ^, 

The access to individual qubits in liquid sample is replaced by simultaneous access to related 
qubits in all molecules of a bulk ensemble. Computers of this type are called bulk-ensemble quantum 
computers. They, in principle, can operate at room temperature. 

For measurements of qubit states the standard NMR technique is used. A principle one-coil 
scheme of experiment is shown in Fig.l. The sample is placed in the constant external magnetic 
field B and the alternating (say, linearly polarized) field b(t), produced by RF voltage V^(t): 

B{t) = B+h{t) = Bk+2bcos{ujt + ip)i, (1) 

where i and k are the unit vectors along the axes x and z. 
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Fig. 1. The principle one-coil scheme of NMR measurement. 



Let the sample be an ensemble of N molecules-microcomputers with L qubits at temperature 
T = 300 K, in the external magnetic fields i? = 1 — lOT. The resonance nuclear spin frequency is 
uja/2it ~ 7ii?/27r < 150 MHz, 71 is gyromagnetic ratio of nuclear spin (71 ~ 7n = 95.8 radMHz/T), 
hcoA/kT < 10-^. 

The output oscillating voltage V{t)) is 

V{t) = QKd^{t)/dt = ^ioQKAdM^{t)/dt, (2) 

where = fioMx{t)dxdz — magnetic flux produced by resonant spins in the coil, Lg = 

IjLo{KA)'^ /V — solenoid inductance of the resonance counter (/io = 47r-10-^T^ cm^/J), Vs — volume 
of the solenoid, K — the number and A — area of coil turns, for resonance oj = uj^ = (LgC)"^/^, 
Q = R/{uJxLs) > 10^ — the quality factor of resonance counter in parallel connected resistance R 
(Fig.l). 

The nuclear spin read-out magnetization (the liquid sample is considered here as continuous 

medium) at optimum resonance conditions is defined by amplitude of RF field b = 1/(71 ^Tj_iT]|i) . 

For transverse T±i and longitudinal T]|j nuclear spin relaxation times ~ 10^ s we have a value 
b ~ 10-3 T and @ 

M,^ax = M,o^r^i/T|,i/2 ^ 7i?i/2 • (N/V,) ■ e(L)/2, (3) 

where N — number of resonant nuclear spins (one in a molecule) in volume Vg- Parameter £{L) is 
the maximum probability of the full nuclear polarization in pseudo-pure state Pj = 1 It may 
be estimated by the difference of equilibrium population between the lowest and highest energy 
states. For nearly homonuclear L spin system it is []•*: 



exp(L^u;A/2A;T) — exp(— LfiwA/2A;T) 2smh{Lfiu:A/2kT) 
^ ~ {exp{hujA/2kT) + exp{-nujA/2kT))i^ ~ 2^ cosh^{hu;A/kT)' 

The other estimations of £{L) are given in |Q. 

In the high temperature limit Huja/ {kT) ^ 1 we have £{L) = L2~^T\uia/ (kT), that is, the signal 
exponentially drops with the number of qubits, but does not for hujA/{kT) » 1 when e(L) = 1 (the 
pure nuclear quantum state). 

^We have used here the two different notation for resonance frequency loa — jaB and loi — jiB. They determinate 
accordingly values e and Mzo in (Q). For homonuclear diamagnetic liquids 7a = 7i- 
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The NMR signal intensity 5 is defined by amplitude 

S = l^maxl = {fio/^)QKA{N/Vshinu;Ae{L), (5) 

The product KA can be expressed also as 

KA = (LsFs//uo)'/' = (Ws/(^oQwa))'/'. (6) 

For the root-mean square noise voltage in the measurement circuit we become 

= VAkTRAiy, (7) 

where Az^ ~ 1 Hz — the amplifier bandwidth. 
So that for signal to noise ratio we obtain 

~ ^J{Q/Vs) • {hujA/kT)Ne{L) • 10"^ (here Vs in cm^). (8) 

For example, for two qubits molecules (L = 2), using, e{L) = hujA/{kT) ~ 10~^, we can make an 
estimate 

(S/N) ~ {Q/V,f/^NlQ-^\ (9) 

Thus, the two qubits liquid ensemble at room temperature for 14 ~ 1 cm^ and Q ~ 10'^ is bounded 
by A'^ > lO^^molecules. 

In the case of paramagnetic liquids one would expect that the number of molecules may be more 
increased with the use of a dynamic polarization (say, Overhauser effect). Assuming for electron 
and nuclear gyromagnetic ratio 7c/7i ~ 10^ we obtain that in the probability £{L) for a L— qubits 
single state the value Tiuja/ (kT) <^ 10^ should be replaced by W^Tiuja/ {kT). Therefore for the same 
value s{L) the allowed number of qubits L approximately will be estimated from 

L2"^ > 10~^ (10) 

whence it follows that L < 12 qubits. 

An additional increasing of read-out NMR signal in paramagnetic liquids may be obtained 
using the ENDOR technique. It is generally believed that for the liquid bull-ensemble quantum 
computers a limiting value is L < 20 — 30 Q . 

However, for realization of large-scale NMR quantum computers a necessary number of qubits 
in quantum register must be L > 10^, that may be attained only at low temperatures in solid-state 
structures. 

The radically new and still unimplemented design of solid-state quantum computer was proposed 
in js], |9| . It involves the formation of an multiple-spin system with access to individual nuclear spins- 
qubits. It was suggested to use an semiconductor MOS structure on a ^^Si spinless substrate, into 
a near-surface layer of which ^^P stable phosphorus isotopes, acting as donors, are implanted in 
form of regular chain. These donors have a nuclear spin 1=1/2 and substitute for silicon atoms at 
the lattice sites, producing shallow impurity states with a large effective Bohr radius. The number 
of donors or the qubit number L in a such quasi-one-dimensional artificial "molecule" may be 
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arbitrary large. It is suggested an electrical control and measurement of qubits states through the 
use of special gate structures. 

The experimental implementation of Kane scheme is undertaken now in Australian Centre for 
Quantum Computer Technology pO[ |. 

There are four essential difficulties in implementing of quantum computer that was suggested 
by B. Kane in gg: 

1) First of all, signal from the spin of an individual atom is very small and it is required of high 
sensitive single- electron measurements. 

2) For initialization nuclear spin states it is required to use of very low nuclear spin temperature 
(~ mK). 

3) In is required of regular donors and gates arrangement with high precision in nanometer 
scale. 

4) It is necessary to suppress the decoherence of quantum states defined by thermal fluctuations 
of gate voltage. 

As an alternative, the variant of an ensemble silicon-based quantum computer was proposed 
by us [0, 12 1 . One would expect that with the ensemble approach, where many independent 
"molecules" of Kane's type work simultaneously, the measurements would be greatly simplified. 
Here we will give a some further development of the scheme considered in [O, 12]. 



1 The silicon structure with regular system of strip gates 

In this case, unlike the structure suggested in , gates A and J form a chain of narrow (/a ~ 10 nm) 
and long strips along which donor atoms at ly distant from each other are placed (Fig. 2). Thus, 
they form a regular structure of the planar silicon topology type. 

The separation between neighbourig donor atoms in Si, as in Kane's scheme, must be < 
20 nm). In this case the interqubit interaction is controlled by gates J. The depth of donor 
occurring d ~ 20 nm. For ly ^ the exchange spin interaction between electrons of donor atoms 
along the strip gates (y— axis) is assumed negligibly small. Then, such a system breaks down 
into an ensemble of near- independent chains-Kane's artificial "molecule, whose electronic spins at 
temperature T < 0, IK are initially fully aligned with the field of several Tesla {'jehB/kT ^ 1). 
The linear qubit density in artificial "molecules" is ~ 50 qubits on micrometer. 

For the initializing of all nuclear spins-qubit quantum states (fully polarized nuclear spins) there 
is a need to attain, for the time being, spin temperature T < 10~^ K. An output signal in this 
system, as in liquids, will be proportional to the number of the "molecules" or donor atoms N 
(components of our ensemble) in the chain along axis y. 
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Fig. 2. The silicon structure (^^Si) with donors atoms '^"'^P and strip gates (Al). 



However for the reahzation of considered structure with donor chains and parallel strip gates, 
as well as of the Kane's scheme, the nanotechnology with resolving of order ~ 1 is also needed. We 
will not discuss here the technology problems in realizing of the suggested silicon-based structure 
(see IC 



2 The states of insulated donor atoms in magnetic fields 

The electron- nuclear spin Hamiltonian for a donor atom '^^P has the form 

H = 7e/iBS - 7i7iBi + ^iS, (11) 

four energy levels of which are given by the well-known Breit-Rabi's formula. For / = 1/2, 5=1/2 
(the z— axis is parallel to B) this formula is written as 



E{F, mp) = - J - lihBmF - (-1)^-^1 + InipX + X^, (12) 

where constant of hyperfine interaction A/{27rh) = 116 MHz 1 13], X = (7e + ^i)tiB/A ~ 'jchB/A, 
F = I± 1/2 = 1,0, and mp = M + m = ±1,0, if F = 1 or mp = 0, if F = (Here M = ±1/2 
and m = ±1/2 are z— projections of electronic and nuclear spins accordingly). The energy levels 
are shown in Fig. 3. For the energy of the ground spin state, F = and mp = 0, hence, we obtain 



F(0, 0) = -A/A - {A/2)Vl + X^. (13) 
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For the next, excited energy state, F = 1, mp = —1 we have 

E{1, -1) = A/4 - (7e - -fi)hB/2. 



(14) 



Thus, the energy difference between the two states of the nuclear spin that interacts with an electron, 
whose state remains unchanged, is described in simple terms (7c 3> 7i for X ~ ^JiB /A S> 1): 



flU) 



A 



A - ^(1, -1) - ^(0, 0) = A/2 + (71 - 7e)?iS/2 + ^ vT+X^ 



71^5 + 



A 



A^ 



2 A-i^hB ' 



hu:-^ = Eil, 1) - E{1, 0) « -j,nB + ^ + 



(15) 



For donor atoms 7e/7i = 1.62 • 10^, -y^ = 176.08 radGHz/T, 71 = 1.137n = 108 radMHz/T. In 
magnetic field B = IT: luX/2tt = 75 MHz, wX/2vr = 41 MHz. 



F 



r(l,l); |1,1) = |l/2,l/2) 



r(l,0); |1,0)= 



(l-c()^^^| l/2,-l/2)+c(^'^|-l/2,l/2) 



r(l,-l); 1 1,-1) = 1-1/2,-1/2) 

£(□,□); !□,□) = 
= (l-d)^^^ I -1/2, l/2)-c(^^^ 1 1/2,-1/2) 

Fig. 3. Energy levels of an individual donor atom in magnetic field. The frequencies 
LOs = ujb — i^A' '^B, f^Ci ^T) = ujb — are in microwave, uj^ - in the RF range of frequencies. The 
transitions with frequencies wb, wd, i^C; i^a allowed, but transitions with frequencies cjs in the 

first approximation are forbidden. 



The state |1,— 1) in M,m basis is equivalent to state \M = —l/2,m = —1/2) and state |0, 0) 
in this basis is the following singlet superposition of states 

|0,0) = (l-a)^/^| - 1/2,1/2) -ai/2|;L/2, -1/2), 

« = ^ fl - -7=7^) (16) 



2 V VTTx^ 

The diagonal matrix elements of nuclear magnetization M^o for two lower energy states will be 
determinated by 

(0,0|M,o|0,0) = (0, 0| 1 0, 0)71 7i= ^=£=71^/2, 

V 1 + 
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(1, -l|M,o|l, -1) = (1, -l|/z|l, -1)71^ = -71^2. (17) 

The populations of the L— qubit lowest and highest energy states for M = —1/2, which give, as 
noted above, the maximum probability of the full nuclear polarisation in pseudo-pure state, are: 

^-^ ^ eM-LfiujX/2kT) 

^' ' (exp(nc^+/2A;r) +exp(-nwX/2fer))-^' 



exp{LnujX/2kT) 
(exp ( hujX /2kT) + expi-hujX /2kT)y 



The maximum nuclear magnetization M^q (the populations of states (1,1) and (1.0) are negligible 
for a;s,WB,u;c » ^X) 



M,o = 71^/2 • {N/V,) 



X eyip{LhujX/2kT) 
Vl + X2 (exp(7lwX/2A:r) + exp(-7ia;^/2A:r))^ " 

exp(-L^a;X/2fcr) 1 



(19) 



{ex.p{hu;X/2kT) + exp{-hu;l/2kT))^ j ' 

For LhujX/2kT <C 1 and X ^ 1 we obtain (compare with (^)) 

M,o ~ 7i?i/2 • (N/Vc) • 2-^L{nujX/kT). (20) 

But for very low temperatures {hujX/2kT ^ 1) we have the full nuclear polarization: M^o ~ 
71^2 • (iV/K) and e(L) = 1. 

3 The gain effect for NMR signal 

Transitions between two lower states are induced by a RF magnetic field, applied at a frequency res- 
onant UJA- The Rabi resonance frequency which is defined by matrix elements of spin interaction 
Hamiltonian with the external RF field h{t) 



Hri= [-fA--flI.)h-b,{t), (21) 

can be found from 

n = -fMX)=2\{0,0\H,i\l,-l)\/h. (22) 
For the amplitude of effective RF field, acting on nuclear spin, 6efj(X) we obtain 

6eff(X) = b (ai/^e/7i + (1 - «)'/') , (23) 

where b — the amplitude of circularly polarized field component. 

The Rabi frequency has the maximum value for X = and rapidly monotonic reduce to value for 
insulated nuclear spin {A ^ 0), 7i6cfr(^) S> 1) = 716. From a rate of quantum operation standpoint 
it is desirable to operate in the relatively weak fields, at which 7e/7i ^ X ~ jc^B/A ^ 1 or 
3.5T > S > 3.9 • 10-3 T. 
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In this case from ( [2^ ) we will obtain 



b,s = (1 + 7])b » b, (24) 

where rj = A/(2^ihB) ^ 1 is the gain factor. On these conditions RF field operate through the 
transverse component of electronic polarization. For magnetic fields B = 1 T we have the value 
6eff = 4.4 • b, and for B = 0.01 T - value bes = 338 • b. 

The gain effect involves an increase of NMR signal and Rabi frequency. In the pulse technique 
it make possible to decrease the length of pulse. This effect was first indicated by Valiev in fl^ . 
The computer operations, owing to this effect, can be done at lower RF fields. At last, it permits 
to reduce the RF field influence on the operation of neighbour semiconductor devices. 

To describe the nuclear dynamics of the studied two low-lying levels system we can write the 
followed Bloch-type equation: 



^!M=^|MxB.,l-^§^-(^^.^. (25) 

where i, j, k are orthogonal unit vectors (Fig.l), 

Beff = (wa/ti) k+26efT cos(a;t)i, (26) 
from where it follows that the value of maximum nuclear read-out magnetization in NMR signal. 



which there have been here for bes{X) = 1/(71 ^Tj_iT]|i), again have value = Mzo/2. That 
is, the read-out NMR signal can not be increased through the gain effect over its maximum value 
Mzo/2, but it can be achieved at more lower RF fields than in the absent of gain effect. 



4 The signal to noise ratio for an ensemble silicon quantum com- 
puter 

For the realization of an ensemble silicon quantum register we propose a variant of planar scheme, 
that contains n ■ p in parallel acting blocks, each have A'o in parallel connected L— qubit Kane's 
linear "molecules". This scheme is schematically depicted at Fig. 4. 

Let the sample be the silicon (^^Si) plate of thick 6 ~ 0.1 cm. For the full number of computers- 
" molecules" in ensemble N = p-N^-n, the volume of sample and also of solenoid isVs ^ 6-lx-ly-L-N 
(the filling factor is assumed to be approximately one). 
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Fig. 4. The scheme of the planar sihcon topology with p ■ n in parallel connected blocks of the 
ensemble L— qubit quantum computers (the connections are not shown here). The broad and 

narrow lines denote the A and J gates. 



The read-out signal from the such ensemble in parallel acting chains, as distinct from liquid 
prototype, for full nuclear polarization or, what is the same, for spin temperatures Tj < 10~^ K 
has no small factor in intensity of the NMR signal of type e(L) = 2~^L ■ Uoj^/ {kT). The NMR 
signal from our sample within of a non-essential factor is the same as from macroscopic sample (see 
Appendix). Therefore, with the expression (^), huJx/{kT) ^ 10^^ (T ^ O.IK) and e = 1 we will 
obtain as an estimation 

(S/N) « ^QhiOA/ikTV,) ■ N ■ 10-9 « ^QN/{5lJyL) ■ IQ-^^. (27) 

It is believed that for superconducting Al gates it is possible Q ~ 10^. The effective volume of 
one "molecule" for l^ = 20 nm, ly = 50 nm, L = 10^, Vs = Sl^lyL = 10~9cm^ we receive that 
the read-out signal in our scheme may be available for standard NMR technique if the number of 
"molecules" in ensemble will be > 10^. High-sensitive devices for individual measurements are 
not needed. 

To estimate the values n, p let us consider the square sample with 50Nqp = 20 • lO^n and 
A'^o = 100. As a result, we receive n ~ 16 and p ~ 63. The area of the structure is ~ 315 x 315 /um^. 
This size is sufficiently small for sample to be housed in the split between the magnet poles of a 
NMR spectrometer. 

For implementation of two-qubit logic operation it is required the controlled by gates J in- 
terqubit indirect interaction with characteristic frequency z^j ~ 100 kHz <^ ujp^l2'K ^ 100 MHz. To 
bring about fault-tolerant quantum computations the relative error for single logic operation must 



be no more then ~ 10~ |15|. Hence it follows that a resolution bound of the NMR spectrometer 
must be of order ~ 100 kHz • 10~^ ~ IHz, that is consistent with the usual requirements. Notice 
that such high precision is needed only for performing the logic quantum operation, but it is not 
needed for read-out measurements. 
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The read-out signal may be some more increased by means of an electron-nuclear double res- 
onance (ENDOR) methods of observing the electron resonance at transition with frequencies 
wb and toc (Fig- 3). 

Consequently, by use of standard NMR and additional of ENDOR techniques the first main 
difficulty of Kane's scheme can be overcame. 



5 The longitudinal relaxation times and long-lived nonequilibrium 
quantum states 

The electron and nuclear longitudinal relaxation times for the allowed transitions in four energy 
levels system of phosphorus doped silicon have been extensively investigated experimentally in ||l( 



see also [17|. For the allowed transitions with frequency ojb and lvc (Fig. 3) electron longitudinal 
relaxation times Tj|B ~ T||c at low temperatures was found to be exceedingly long. They are of 
order one hour at T = 1.25 K, B ~ 0, 3T, are independent of phosphorus concentration below 
C ^ 10^^ cm~'^ and are approximately inversely proportional to lattice temperature T. But they 
very fast drops to ~ 1 s already at C ~ 10^''' cm~^. For the structure on Fig. 2 d = 20 nm, l^ = 20 nm, 
Ij = 50 nm we obtain the average distance between atoms ^^P in silicon Zav = (20-20-50)^/^ = 27nm 
and their associated concentration C = (/av)~^ = 5.10^^ cm. For this value of C according to [^] 
we must have Tj|B 100 s. 

The further increasing the electron relaxation times may be achieved by decreasing of lattice 
temperature. It is possible to extrapolate electron relaxation time to T ~ 0.1 K resulting in a value 
7||B ~ 1000s. 

The relaxation time for transition with frequency loj), which involves a simultaneous electron- 
nuclear spin flip-flop, at T = 1.25 K, C ~ 10^^ cm~^ and i? ~ 0, 3 T was T^^d ~ 30 hours » T]|b, T]|c. 

The nuclear longitudinal relaxation time T||j exceed 10 hours |jl^. Notice, that for full electronic 
polarization, occurring at T < 0, 1 K, this time does not depend on the interaction with neighbour 
electron spins [1T7|| . 

The dipole-dipole interaction of separated qubits in our chains with polarized electron spins 
of neighbours phosphorus atoms is divided into secular and nonsecular parts. The main part 
of it is the secular part. It lead up to constant shift of resonance nuclear frequencies of order 
/^o7N7B^/(4'/r/^) ~ lOHz, where 7b = 2^b/^ = 176radGHz/T, l^ ~ 20 nm. 

To suppress the influence of the unwanted secular part of dipole-dipole spin interaction in a 
chain, may be applied the external magnetic field, B, at an "magic" angle, 9 =arccos(l/\/3) = 
54.74*^ to direction of the chain. 

The extremely long relaxation times of the electron and nuclear spins imply that the required 
initializing of nuclear quantum states (full nuclear nonequilibrium polarizations) can be attained 
by deep cooling of short duration only nuclear spin system to T < 1 mK without deep cooling of 
lattice. There is the possibility to reach it at indirect cooling of nuclear spin system by means of 
dynamic nuclear spin polarization techniques [17|. 



6 The cooling of nuclear spin system and initialization nuclear 
states by means of dynamic polarization 

One such method of dynamic nuclear spin polarization for donor atoms is based on the saturation 
by the microwave pumping of the forbidden transition (frequency ws in Fig. 3), that is designated 
as the solid state effect ||6|, |l^ . 
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Let us consider this effect. Tlie polarizations of electrons Ps = '2,{Sz) and of nuclei F\ = 
can be expressed as 



Ps = 1) + 0) - -1) - p(0, 0), 



(28) 



Pi = p(l, 1) + p(0, 0) - p(l, 0) - p(l, -1) 



where p(P, Mp) are the populations of states P, Mp (Fig-3). They fill also the requirement 



dp(0, Q)/dt = (p(l, 1) - p(0, 0)rB)/2||B + (p(l, 0) - p(0, 0)rD)/T||D+ 
+ (;,(!, -l)-K0,0)r+)/7||A, 
dp(l, -1)1 dt = (p(l, 0) - p{l, -l)rc)/2||c + (p(l, 1) - p(l, -1)) • 

+ (p(0,0)r+-p(l,-l))/2||A, 
dp(l, 0)M = (p(l, -l)rc - p(l, 0))/7||c + p(0, 0)rD - p(l, 0))/2]|d + 

+ (p(l,l)-p(l,0)rX)/2||A 
dp(l, 1)M = (p(0, 0)rB - p(l, 1))/7||B + (P(l, -1) - 1)) • We+ 



where parameters rB,c,D,A = exp(— ^u;B,c,D,A/fcP) are ratio of rates for a up and down thermal 
transitions. For values hcOB,c,r)/kT » 1, hu^/kT <^ 1 (T < 0, IK) there are the thermal electron 
Pso = —1 and nuclear Po = Tiuj\/kT <^ 1 polarizations. 

Let us propose next that rate Wc of the induced forbidden electron transitions |1, 1) =^ |1, — 1) 
and electron longitudinal relaxation times are satisfied the conditions: < T]|b ~ T]|c,T]|a ^ 

Tjp < Tj|g(rj|s is the longitudinal relaxation time of electron spins for forbidden transition). Here- 
after we shall write 



p(l,l)+p(l,0)+p(l,-l)+p(0,0) = L 



(29) 



Let us write the rate equations for the populations: 



+ (p(l,0)rX-p(l,l))/T||A 



(30) 



dp{0,(})/dt 



p(l, 1)/7||B + ip{l, -1) - p{0, 0))/7||A, 



dp{l,-l)/dt 



p(l, 0)/7||B + ((p(l, 1) - p{l, -1)) • We+ 



+ (p(0,0)-p(l,-l))/2||A, 
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dp{l,0)/dt = 



-p{i, o)/r|,B + {p{i, 1) - p{i, o))/r|,A+ 



dp{l, l)/dt = -p{l, 1)/2|,B + -1) - p(l, 1)) • VFe, 

+ (p(l,0)-p(l,l))/2||A (31) 
With equations (p8|),(p9|), (l3l|) we can obtain the rate equations for and Py. 

dPs/dt = -{Ps + Pi) ■ We - {Ps + 1)/2||B, 

dPi/dt = -{Ps + Pi) • We - PimA. (32) 

The steady-state saturation condition {We ^ 1/T\\p^) of the transition |1,1) =^ ll^— 1) gives rise 
to the equahzation of the populations p(l, 1) = p(l, —1) and respectively to the full nuclear spin 
polarization 

Pi = -Ps=p{0,0) = l. (33) 

It is obvious that this state is equivalent to state of cooling up spin temperature Ti < hujx/k ~ 
10~^ K nuclear spins. 

Let us estimate now the needed microwave power for saturation. The rate of external microwave 
field that induce forbidden electron transitions with frequency (for X k, ^e^B/A < 1) is 

We < (7c6mw)' • (A^)-'/', (34) 

where (Awg)^/^ ~ ^T^g"^ — nonhomogeneous broadened resonance line widths for the saturated 
electron transition, is amplitude of microwave field. The quality factor Qe of microwave cavity 
is 

Qe-u;sbi^-Vr/{2fioP), (35) 
where — volume of microwave resonator, P — dissipating power. The saturation condition is 

We » l/2]|s or {^ebm^fTlsTiis » 1, (36) 
The power dissipated in cavity is 

P > LOsVrWe{A^)'/^/{2floQeJ^). (37) 

For example, taking We ~ W^s''^, ws ~ 100 radGHz, ~ Icm^, Qe ~ 1000 and (A^)^/^ ^ 
10^ s~^ jl^ as a rough estimate we obtain P > ImW. Notice that this power is applied only 
during the saturation process over of time >W-^^l ms in the act of initializing of qubit states. 

Hence, the initialization of nuclear states may be obtained by use of ENDOR technique at the 
lattice temperature of order 0,1 K and by this means the the second difficulty of Kane's scheme can 
be overcame. 
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7 A silicon-based ensemble NMR quantum computer of cellular- 
automaton type 

For implementation of ensemble silicon quantum computer, operated on cellular-automaton princi- 
ple, it may be usable the previously considered system of long chains of donor atoms ^^P disposed 
in silicon, but free of the A and J gates. 

If exchange interaction constant (it is here positive) for localized electronic spins of ^^P along 
the chain is more than Zeeman energy J{1) ^ jb^B and the temperature is well below the critical 
temperature for electron ordering (Neel temperature) T < Tns ~ J{l)/k {k = 1, 38-10^^^ J/K — the 
Boltzmann constant), the one-dimensional antiferromagnetically ordered ground state of electronic 
spins is produced. We can estimate the exchange interaction constant as J > 7N?i-B ~ 6.5 • 10~^^ J, 
that correspond to the distance between donors 1^, ~ 20 nm |Q, that is the electron spin critical 
temperature will be Tns ~ 4K. 

Due to hyperfine interaction nuclear spins will be oriented according to the electronic spin 
direction in the resultant field and will form array with the alternating orientation of nuclear spins. 
At magnetic fields B < A/^^h ~ 3.5 T and at temperatures T ~ 10~^ K the nuclear spins ^^P will 
form a periodic ground state array of ABAB. . .type: titi . . ., where | marks the ground state of 
nuclear spin in an A-site and [ — the ground state of nuclear spin in a B-site with almost 100% 
opposite orientation (u;a,b/^T > 1). That is the distinct nuclear spins are in initializing ground 
state. 

Notice that the using of dynamic methods, makes possible the high orientation of nuclear spins 
also at larger temperatures and, according to the above section 5, this state will be the long-lived 
nonequilibrium nuclear spin state. 

The nuclear resonant frequencies wa,b of neighbor nuclear spins are different for each of the 
magnetic one-dimensional subarray A and B in the chain as they depend on the states of neighboring 
spins. We will take it here in the simple form |12, |l9| ]: 



WA.B = w(m< -I- m>) PS l^ihB ± A/2 - /n("i< + m>)|/?i, (38) 

where In ~ {A^ / J) is constant of two neighbor nuclear indirect spin-spin interaction due to hy- 
perfine interaction, m< and m> are the magnetic quantum numbers for the left and right spin 
subarray. The nonsecular part of nuclear-nuclear interaction is neglected here taking into account 
that 7n^-B, a/2 ^ I^. The difference between nuclear resonant frequencies for the distinct neigh- 
boring spins orientations is Aa'i/27r ~ In/2TTfi ~ {A^ / J) /2iTh ~ 0.5 MHz, whereas the resonant 
nuclear frequencies are a;A,B/27r ~ A/iirh ~ 120 MHz. 

For the organization of logic operations let us use the addressing to spin states, analogously to a 



scheme put forward in [20|. Each nuclear spin in A-site of this scheme has two internal eigenstates 
- ground | t) and excited | JJ.) and in B-site, accordingly, - | [) and | ft). 

We take into account that the life time of excited states (the longitudinal nuclear spin relaxation 
time Tjj) at low temperatures is very long. Each logic qubit of quantum information in this state 
will be encoded here, similar to |2^, by the four states of physical spin-qubits (individual nuclear 
spin states): the logical qubit basis state "0" will be encoded by unit | -Ij-ftti)) while the state "1" — 
by I tiJJ-'fr)- It is important here that the resonant frequencies of nuclear spins depend on neighbor 
spins states. Both logical states have two excited spin states and the zero projection of total nuclear 
spin. 

Notice that a random inversion of only one spin will result in degradation of the qubit state. 
But to form the rough error, for example, of "0" ^ "1" type in the coding of stored quantum 
information it is essential to invert simultaneously four spins. Therefore, it may be concluded that 
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the considered way of qubit coding ensures a better fault-tolerance with respect to this type of 
errors. 

As the ports (noted on Fig. 5 by x) for input and output of the information in the array of 
ground states spin-qubits can be dopant nuclei D at the certain place of the array with distinct 
resonant frequency, defects or local gates that modify the resonant frequency of the nearest nuclear 
spin in the array. 



Electron spins 5 

Nuclear spins I tjtjxl^ 

port 



Fig. 5. Scheme of electron and nuclear spins ordering 



Starting from the perfectly initialized states inputting the information can be performed by 
setting the dopant D-spin to a desired state by means of RF-pulse at its resonant frequency. The 
nuclear spin state of spin nearest to the dopant spin is created by SWAP operation. After the 
required information is loaded, D-spin is reset to the ground state. Upon completion of computation, 
the state of any spin can be measured by moving it to the A-site nearest to D, then swapping A=^ D 
and finally measuring the state of D-spin. 

For the implementation of quantum operations on logic qubits we will also introduce, as in (p!^), 
one auxiliary control unit (CU), which is represented here by six physical spin states in the pattern 
ft-IJ-iTft-ll-- The CU exists only in one place along the array and is separated from logical qubits 
by odd number of spacer spins. The applying corresponding SWAP sequence of pulses CU leads 
to interaction of CU with one and two logical qubits and performing on they one- and two-qubit 



quantum operations (for more details, see [|12|, |19[). 



In the case of a large enough ensemble of in parallel acting chains the states may be measured 
by NMR methods. For the increasing of logic qubit number in sing le "molecule" to > 10^, that falls 
on one port, it can be used two- and three-dimensional structures with antiferromagnetic chess-type 
ordering of electron spins. The corresponding ordering will have and for nuclear spins (Fig. 6). 



port 

1 t 1 M t 1 t 1 t 1 ^ 



t i t i t i t i t i t - 



t i t i t i t i t i t 

Fig. 6. Scheme of two dimensional ordering of initialized nuclear spins. 
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If, for example, the number of spins-qubits, falls on one port in linear chain is, say, L ~ 30, so in 
two-dimension case their number will be L = 900. Ensemble that is composed of iV ~ 10^ in parallel 
acting such plane "artificial molecules" permit to provide input and output of information through 
the standard NMR techniques. The using of more sensitive ENDOR techniques has particular 
meaning when it may be combined with the techniques of dynamic polarization (solid state effect). 

The advantage of ensemble quantum cellular automaton in the comparison with the above 
considered ensemble variant with strip gates is as follows: 

a) The system of the control gates is absent, what essentially simplifies the producing of com- 
puter structure and eliminates one of the important source of decoherence. 

b) The coding of logic qubits into four physical qubits gives a higher degree of fault-tolerance 
in logic operations. 

It follows that the going to ensemble quantum cellular automaton permits to overcome the third 
and fourth difficulty of Kane's scheme. 

The main disadvantage of cellular automaton scheme is the relative complexity of logic operation 
performing. 

Notice here that there is also an another possibility of the ensemble NMR implementation, 
which does not have the gate system. The selectivity of nuclear resonance frequencies for individual 
qubits in the chain ensemble can be achieved, instead of A— gates voltage, by applying the external 
magnetic field gradients along axis x for neighbour qubits separated by ~ 20 nm of order dB^/dx ~ 
IT/cm (that is feasible now), which produce a ~ 100 Hz resonance frequency difference. 

8 Decoherence of the qubit quantum states in an ensemble regis- 
ter 

The interaction of nuclear spins-qubits with polarized random distributed paramagnetic centre 
(atoms and defects) in substrate causes a spread of resonance frequencies for distinct qubits, deter- 
minated by the dipole-dipole electron-nuclear spin interaction. For low temperatures (T < 0, 1 K) 
it is the main mechanism of nonhomogeneous broadening of nuclear resonance line widths (Awj )^/^. 
Therefore for its elimination it is required to use the cleaning from such centres. 

To estimate the allowed concentration of paramagnetic centres in silicon substrate Cs we will 
take the expression 

(A^)^/^ ~ //o/(47r) • 7N7B?iCs > l/T^i » 1/T\\i, (39) 

where T±i — the homogeneous nuclear spin transverse relaxation time, which correspond to deco- 
herence time td of nuclear quantum states. This time to be sufficiently large (~ 1 s) as compared 
with the RE control pulses length (~ 10"^ s) in order one can be done in a time of order td a 
necessary for large-scaled computing number of logic operation (> 10^) ]T5[ . 

Taking into account that frequency of interqubit indirect interaction z^j w 100 kHz (see above), 
for the coherence length of quantum register state we obtain A ~ 27rAj/i/x • « 6.3 • 10^ • l^- 

With Td = 1 we obtain for the allowed concentration of paramagnetic centers the estimate 
Cs < 10^^ cm~^, that is almost usual requirement for a degree of cleaning from unwanted impurities 
in silicon technology. 

The another mechanism of decoherence of qubit states is the dipole-dipole interaction with 
nuclear spins of impurity atoms of nuclear spin containing elements. 

The allowed concentration Cn of impurity with nuclear spins in silicon substrate can be esti- 
mated from the following expression 
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(40) 



Taking again td ~ Is, we will obtain Cn < 2 • 10^^ cm ^. For the number of silicon atom 5.0 • 
10^^ cm~^ the allowed relative concentration of the impurity nuclear spin, contained isotope ■^^Si, 
must be Cn% < (2 • 10^^ • 100)/(5.0 • lO^^) = 0.4 • 10"^%. For comparison, the natural abundance 
of isotope ^^Si in natural silicon is 4.7%. At present the realized degree of cleaning ^^Si is 99,98%, 
which is yet not sufficient to our purposes. 

The allowed concentration for another impurity elements having the natural abundance of 
isotopes with nuclear spins P < 100% is determined from 

Cn% < IQ-^WO/P. (41) 

In the case of ensemble approach for implementing of logic quantum operations, for read-out of 
information and for suppressing of decoherence effects, caused by secular part of dipole-dipole 
interaction qubits with randomly distributed paramagnetic centres, may be explored the well known 
many-pulsed technique of high resolution NMR spectroscopy in solids [^]. 

As source of errors in quantum computer operation there is also generation processes unwanted 
states for neighbour initialized qubits as a result of nonresonance effects. This is a further mech- 
anism of the decoherence of quantum states. The probability of such errors can suppress to the 



required value (< 10 ) by choosing of Rabi frequency value for control RF-pulses |22]. For one 



example, if resonance frequencies of nonresonant qubits differ from the frequency of resonant qubit 
by Auj, so the Rabi frequency for vr— pulse should be = \Auj\/\/Ak'^ — 1, where fc = 1, 2, 3, . . . 

Notice that it was also suggested another, as compared to considered above, scheme of the 
four-spin encoding two logical qubits, which are represented by the two zero-total spin states of 
four spins, generated by the pairs respectively of the singlet and triplet states [^]. This scheme 
is based on the so named collective decoherence conditions. This conditions can be attained in 
coupled spins at very low temperatures, when the longest wavelength magnon modes are quenched 

MM- 

The chief disadvantage of the variant with strip gates as and Kane's scheme is associated with 
availability of the gates. The thermal voltage fluctuations on gates are additional mechanism of 
quantum state decoherence. They are essentially suppressed here if both the Al gate electrodes 
and the entire resonance circuit are superconductors at the work temperature T <Tq {Tq = 5.8 K 
for films). 



Conclusion 

1) The line of the large-scale ensemble NMR quantum computer development have certain advan- 
tage over the line, based on Kane's scheme. It consists in the possibility of employment of the 
standard NMR technique for the measurement of quantum states at output of computer, like in 
the liquid prototype. 

2) For the initialization of nuclear spin states at temperature T ~ 0, 1 K methods of dynamic 
polarization of solid state effect type may be proposed. 

3) Analysis of proposed planar structure of ensemble silicon computer shows the possibility of 
realization of large-scale NMR quantum computer for ensemble component number N ~ 10^. 

4) The implementation of cellular automaton principle permits to abandon the realization reg- 
ular nanostructure in form of gate chains. 
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5) Analysis of different feasible ways for to obtain a decoherenec times large enough shows that 
the values, needed to perform required for large-scaled computations number of quantum logic 
operations ~ 10^, can be achieved. 

6) One would expect, that in the future the combined variants of solid state quantum computers 
will be proposed. They can, say, exploit the structures that are contained nuclear spins and quantum 
dots with electron spins and the combined methods access to qubits, like ENDOR method. 



Appendix 

Signal NMR for discrete ensemble of nucleeir spins 

Let us consider here the sample that involves N = nL ■ pNo nuclear spins arranged in the plane 
z = of the silicon plate at regular intervals along the strips (Fig. 4). The spins in chains under 
strip at resonance in the each block are oriented along x axis (solenoid axis) and separated at 
intervals of L. The read out NMR signal is 

K /•^/2 f 

\Vjnax\ = QuA-rp / I / B^rna:K{x,y,z)dydz\dx, (42) 
X J-X/2 J A 

where X = nL ^ L is roughly the length of solenoid and 



-Bxmax(a;,?/,z) = (43) 
^ -2{x - Ln,f + {y- l,p,f + 

is the peak magnetic field produced by resonant spins in solenoid. For simplicity it is suggested 
that n,p,NQ,L 1 are the even numbers. We assume that the area of coil turns is A = D ■ 6 

{D = ly-pNo, S-^D). 

For summation over nj and pi we have used the Poisson's summation formula, namely. 



pNo/2 



E /(^yfi) = ^E / meMi^^</ly)d^, (45) 



Pi=-pNo/2 v=-oo' 

by omitting the oscillated terms with u ^Q: 



|Vmax| = jJioQKiOA-^^ / dx— / dr]- 

OTT A J-X/2 X J-X/2 

D J-D/2J-D/2[{x-r]f + {y-£,Y + z^f/^ 

-)2 ^ ^2 



Taking into account D z , upon integrating (A.4) over y, ^, we obtain 



StT D X J-X/2 X J-X/2 
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By integrating now over z and preserving only the major logarithmically increasing for 2\x — rj\/5 =^ 
term, we obtain 

\Vrae,yi\'^ ^iGQKAuJA-, TTr^ / "^"F / log ■ 



47r AD X J-x/2 X J-x/2 \x — ri\ 

IV X X 

= {f,o/^)QKAiOATT7ih ■ — log (48) 
Vg itL> dy/e 

where K = AX, N = npNo, e = 2,718. . . . 

We see that expression ( ^8[ ) is distinguished from (Q) by non-essential factor ^ log , that 
is of order several ones. 
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